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LONG-TERM STORAGE CAPACITY 
OF RESERVOIRS 


By H. E. HURST? 


SYNOPSIS 
A solution of the problem of determining the reservoir storage required on a 
given stream, to guarantee a given draft, is presented in this paper. For ex- 
ample, if a long-time record of annual total discharges from the stream is avail- 
able, the storage required to yield the average flow, each year, is obtained by 
computing the cumulative sums of the departures of the annual totals from 
the mean annual total discharge. The range from the maximum to the mini- 
mum of these cumulative totals is taken as the required storage. 
The results of the investigations given here were applied in computing the 
storage required in the Great Lakes of the Nile River Basin, the length of time 
chosen as a basis of the estimates of storage being 100 years. These reservoirs, 
in combination with the other projects described, would enable the Nile to be 
developed for irrigation to the fullest possible extent. Finally, a relation has 
_ been derived between the storage capacity on the main stream and the amounts 
' required on the tributaries to produce its equivalent. 
It is thought that the general theory may have other applications than the 
- design of reservoirs for the storage of water. 


1. INTRODUCTION 


A number of projects has been investigated which will enable Egypt and 
the Sudan to develop irrigation from the Nile River to its fullest extent.?* 
The broad plan involves storage of water from good years for use in bad ones, 
and necessitates reservoirs of sufficient capacity to meet the shortages that 
might occur during a century. Much research had been done on the capacity 


Nore.—Written comments are invited for publication; the last discussion should be submitted by 
September 1, 1950. 

1 Scientific Consultant to the Ministry of Public Works, Egypt; formerly Director-Gen., Physical 
Dept., Cairo, Egypt. 

2'‘The Nile Basin,” Vol. VII, by H. E. Hurst, R. P. Black, and Y. M. Simaika, Physical Department 
Paper No. 51, Ministry of Public Works, Cairo Govt. Press, Cairo, Egypt, 1946. 

3 Civil Engineering and Public Works Review, September, 1948. 
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of such reservoirs and the results were used to determine the size of the reservoirs 
required in the Great Lakes of the Nile Basin for what was called “century 
storage.” . 

Some preliminary results of investigations were released in 1938* which 
were amplified in 1946;? but in neither case was there an adequate statement of 
the theoretical investigation. The latter has been reserved for this paper. 

Storage problems were treated first by the late Allen Hazen,’ M. ASCE, in 
1914, when he studied the discharge of thirteen American rivers, dealing with 
them mainly in a graphical manner, on probability paper. The longest record 
was one of 45 years, but longer sets of observations were made by combining 
records from several rivers. The objection to this procedure and the uncer- 
tainty in extending the normal Gaussian frequency distribution to the future 
were recognized by Mr. Hazen, but he thought it gave the best results then 


~ obtainable. 


The work was extended by the late Charles E. Sudler,* M. ASCE, by similar 
graphical methods. Much of Mr. Sudler’s work was based on artificial records 
in which information taken from a short record of a stream was extended to a 
period of 1,000 years by writing, say, 50 annual runoff values on cards and 
then shuffling and drawing a card from these 1,000 times. Again the actual 
data were necessarily confined to short periods since stream measurements in the 
United States had rarely been made for as long as 50 years. Thus, the results 
are mainly based on a few short-period records of natural river flow extended to 
long periods by random repetition. It is clear that the process can give very 
little information about what may happen over long periods and may be com- 
pletely misleading since it excludes values higher or lower than those which 
have already occurred. This was recognized by some of those who contributed 
to the discussion on Mr. Sudler’s paper. Moreover, the graphical methods 
employed do not lead to any concise and easily understandable presentation 
of results. 

In the present investigation the methods of the theory of probability are ; 
applied, but the data are long-period records of natural events which are shown 
to have certain similarities, although they may be records of river discharge, 
rainfall, temperature, annual growth rings of trees, or annual deposits of clay 
in lakes. The average results can be expressed quite simply by two equations 
which apply to the storage of other substances as well as water. It is shown | 
that, in regard to storage, records covering such short periods as 30 years or : 

40 years may be very misleading, but their usefulness can be extended by 
analyzing a large number of phenomena which have been recorded for periods — 
as long as a century or more. io an 

The problem of what storage is required on a stream to give a certain 
minimum discharge was first investigated in the case of Lake Albert (in the 
Belgian Congo and the Uganda Protectorate) to determine the size of an over- _ 
year storage reservoir which would be required in order to equalize the outflow 


46 MY +f ” 
4**The Nile Basin,’’ Vol. V, by H. E. Hurst and P. Phillips, Cairo Govt. Press, Cairo, Egypt, 1938, p. 86. _ 


5 “Storage to Be Provided in Impounding R: i ici i 
Transactions, ASCE, Vol. LX XVII, Desaettien, 1. bis Water. Sopris. Se aaa ) 


i ae Required for the Regulation of Stream Flow,” by Charles E. Sudler, tbid., Vol. 91, 1927, 


- ? 
ee 


RESERVOIR CAPACITY 5 


over a number of years. It leads to a statistic R, which is the range from maxi- 
mum to minimum of the curve obtained by plotting the cumulative totals of 
the departures from the mean of the annual discharges taken in order. Such 
} a curve of cumulative totals is sometimes called a mass curve. Values of R 
have been computed for long series of observations of rainfall, temperature, 
and other natural phenomena, and these have been filed for reference in the 
Engineering Societies Library.® 

Nomenclature——The letter symbols introduced in this paper are defined 
‘where they first appear, in the text, or by illustration, and are arranged alpha- 
betically in the Appendix, for convenience of reference. A few specialized 
terms used in the paper are defined as follows: 


Century Storage——The computed storage capacity of a reservoir, as re- 
quired in the Great Lakes of the Nile Basin, necessary to meet the shortages 
that may occur in 100 years (Section 1). 

Continued Sum.—Successive cumulative totals (Section 1). 

Cubic Meter.—English units, 0.0008107 acre-ft. 

Equatorial Lakes——Lake Victoria (in Uganda, Tanganyika, and Kenya), 
Lake Albert (in the Belgian Congo and Uganda), and other bodies of impounded 
water in the southern Nile Basin, at the equator; also referred to as the ‘‘Great 
African Lakes,” or simply the ‘‘Great Lakes” (Section 2). . 

Gaussian Frequency Distribution—The normal frequency curve first advo- 
eated by the German mathematician, Karl Freidrich Gauss, who maintained 
that all skew distributions tended to become normal as the number of events 
was increased and that skewness was an evidence of inadequate data. 

Meter.—English units, 3.28083 ft. 

Milliard—One thousand million things—that is, 10°. 

Varves.—Annual layers of silt as deposited in a lake or other body of still 
water. As individual layers differ in thickness and character, a succession of 
such layers forms a characteristic group which can be identified as of con- 
temporaneous deposition in whatever deposit it may be found. It is thus 
possible, by combining different sections, to measure the time involved in the 
deposition of the entire group of sediments and to construct a time scale in a 
manner similar to that employed in the study of annual rings in trees (Webster’s 
_ New International Dictionary of the English Language). 


2. STATEMENT OF THE PROBLEM 


Given a reservoir of large capacity, how can its outflow be regulated eco- 
nomically to meet the deficiencies of low years, and what is the capacity required 
to guarantee a certain minimum discharge in these years? These are the 
basic questions. 

As a concrete example to illustrate the methods, consider the case of Lake 
Albert. If the past annual discharges from the lake are given, the size of 
reservoir required to maintain the maximum possible steady discharge during 
the period can be determined. This is the mean discharge for the period, and 
the capacity of the reservoir is obtained by adding, algebraically, the de- 
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partures from the mean to form a series of accumulated totals. The difference 
between the highest and the lowest of these accumulated totals—that is, the 
range R—is (a) the maximum accumulated storage when there is never a 
deficit, (b) the maximum accumulated deficit when there is never any storage, , 
or (c) their sum when there is both storage and deficit. This range R is the 
storage required to maintain the average discharge. 

Increased losses due to storage are disregarded because, unless they are 
small, the site is not suitable for over-year storage. On a reservoir like the» 


Annual out flow. 


80 Continued sum of 
departures. 


AA, Mean annual out flow. 
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Milhiards of cube metres. 


1903 1910 1920 1930 1940 1946 


Fie. 1—Annvat OvuTFLow FROM LAKE ALBERT IN THE BELGIAN ConGco AND UGANDA 
PROTECTORATE (1904 Tro 1946) anp THE CuRVE oF ConTINUED SuMs oF 
Departures (OnE Miuirarp Cusic Meters = 810,700 Acre-Ft) 


Aswan Reservoir in Egypt there is heavy evaporation uncompensated by any 
rainfall, and the surface area of the reservoir is several times that of the natural 
river. Therefore, there is an appreciable extra loss as a result of storage, so 
that, even if there were the capacity at Aswan, it would be unsuitable for over- 
year storage. In the case of the Equatorial Lakes their size makes it possible 
to store a great volume of water with only a small proportional increase of 
surface area. Moreover, evaporation and rainfall are nearly equal; so in- 
creased evaporation loss due to increased surface is compensated by increased 


gains from rainfall. 
% 
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Fig. 1 shows the annual total discharge from Lake Albert based on recorded 
lake levels and gage-discharge curves. The curve of accumulated departures 
from a base near the mean is also shown. The actual values and the method 
‘of computation are shown by reference to Table 1. 


The discharges, Q, in Col. 2, are 
usually given to one tenth of a TABLE 1.—DiscHarce FROM Laker 


milliard cubic meters but the last ALBERT (BrLcian Conco—UGanpa 


digit is not certain. For this com- PROTECTORATE) IN MILLIARD CUBIC 
putation the discharges are written Meters (1 Miuuarp Cusic 
to the nearest milliard and this re- Mertmrs = 810,700 Acre-Fr) 


duces most of the departures, d, to 


a single digit, which is enough for BA Soha TDN IS weal weeeeed Plea 
purposes of computation. The omis- GP POs Ce 1G = Oe ke) 
sion of the last digit makes no sig- ee a0 H paras 
nificant difference to the mean or the 1906 34 10 100 |. 28 
ae 1907 33 9 80 | 37 
standard deviation, o, but lessens 1908 26 2 0} 39 
considerably the labor of computa- : ae 
tion. The computation is arranged —jgit x 3 . Oi ieae 
b : : 1912 19 5 20 | 39 
to be self-checking as far as possible. 1912 aa a at mets 
The mean of the series is first tte a : 10 a 
computed, and then a base is chosen 1916 2 3 10 | 35 
: 1917 47 23 530 | 58 
near the mean. In this case the 1918 48 24 580 | 82 
mean is 1,040/43, or 24.19, and 24 is - ae ee ae! 
selected as a base from which to 1931 1 bream Rees 
. cf it 11 | 120] 68 
compute departures, being the near eee re Ay he cob: (eles 
est whole number to the mean, and ree ae g aE 
so keeping all the departures to one 1925 19 ; 5 20 39 
or two digits. 1928 21 3 10 | 37 
After writing the departures the z 
positive and negative values of d ee aa 2 : slices 
(Cols. 3 and 4) are added separately. ee 28 : Sole as 
_ In this case the sum of the positive is iets 23 : on A 
122 and of the negative, 114. The 1936 20 4 20 | 31 
; ; 1937 24 0 O.lcsat 
difference is 8 and 8/N should be 1938 26 2 ‘ | 33 
the difference between the mean and 
, : : 2 4 
thebase. This latter difference gives soy 2 an Parte NW 
a check on the correctness of the'de- $43 ae etree 
; ; 1944 18 6 40 | 25 
partures, and also on their continued Tae a alles ie 
sum 2d, (Col. 6). If the departures 1946 | 16 8 60 8 
were reckoned from the mean their z 1,040 | +122 | —114 | 2,310 


final algebraic sum would be zero; Mean dato | 048 
otherwise, it is N (mean minus 
base), or 43 (24.19 — 24), and, of 
course, this is also the algebraic sum of the totals of the positive and negative 
departures, d, thus affording a check on the continued sums. With a long 
series of observations it is advisable to compute the totals of the columns, 2, 
on each page to check the work piece by piece. 
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In computing the standard deviation, ¢, the mean square deviation must it 
be reduced by (mean minus base)*, or by (24.19 — 24)? = 0.19? = 0.04, to give : 
= 2a — 0.04; and o = 7.33. This is the only part of the work which is 
not self-checking and must be checked by repetition. The last digits of the 
square are omitted, which means in this case that the squares are given to the 
nearest ten which is near enough provided there are 50 or more observations, 
The effect of the approximation on o is about 1% in this instance. 

The departures are reckoned from a base which is not the mean, and, if the 
continued sums of these are plotted, a simple graphical method will yield the 
‘continued sum of departures from the mean or any other base. In the present 
case the mean is 24.2 and the base, 24. Every departure from the base is 
therefore 0.2 greater than the corresponding departure from the mean, so that 
the sum of the first n departures is n X 0.2 greater than the sum of the first n 
departures from the mean. The sum of all the departures from the base is 
N X 0.2. Thus, the continued sum of any number of departures from the 
mean can be found by joining the point representing the last continued sum to 
the origin and taking this line as a new axis. This is a useful construction for 
finding the continued sums from any base, when the sums referring to some 
base have been plotted. 

In Fig. 1 the continued sums of departures from base 24 are plotted and 
the new axis OX giving continued sums of departures from the mean is drawn. 
The range, R, of the summation curve of departures from the mean is PM = 84 
milliards. 

Examples of other uses of the graphic method of obtaining continued sums_ 
are given, subsequently, in Figs. 3 and 6. 

In the case of Lake Albert the accumulated sum of departures rose to a 

‘maximum in 1919 (point P, Fig. 1) and then fell to zero by 1946 (point X), so 
that the next period would begin with no stored water. Had the operation 
_ been begun in 1920 there would have been a deficit to start with. This fact 
shows the necessity of making a reservoir some time before it will be required 
so that it can be partly or wholly filled before its storage is needed. It is clear 
from Fig. 1 that the storage required to maintain the mean discharge over the | 
period is the range from maximum to minimum of the accumulated departures. 
It is also clear that, if the capacity of the reservoir is just enough to make up | 
the maximum accumulated deficit, when this has been done, the reservoir will 
be empty and there will be no margin of stored water for the future. Moreover, © 
the engineer is regulating the flow after the event, with full knowledge of what 
has happened. The problem of the required capacity is not a simple one and 
cannot be solved by adding departures from the mean, considering only Lake 
Albert, but must be dealt with on a wider basis. ‘The discharge of the main © 
Nile River at Aswan may be taken as a further illustration of the difficulties. 
Table 2, for example, is based on three widely differing assumptions: First, , 
that a computation was made in 1909, using 38 years of record; second, that a 
computation was made in 1946, using the succeeding 36 years of record; and, © 
third, that a computation was made in 1946, using the total period of 74 years _ 
of record. The difference between the results of these computations shows. S 
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how misleading the simple process can be, both as to the size of the reservoir 
‘needed and as to the benefit to be expected fromit. Thus, it is clear that more 
‘extended investigation is required, and this can only be made by the examina- 
“tion of many sets of river discharges extending over long periods, or of other 
data with similar characteristics. At the same time some promise of a solution 
is possible by exploring theory. 


TABLE 2.—Comparison OF VARYING RESULTS, COMPUTATIONS FOR 
THE Main NILE RIVER 


Description 1871-1908 1909-1945 1871-1945 
Average annual discharge, Q, milliards m3............. 103 83 b 93 
Range of accumulated departures, R, milliards m’...... 201 83 476 


LEE 


The first part of this paper attempts to find a theoretical relation among the 
range, R, of the summation curve of departures from the arithmetic mean, the 
variability of the phenomenon expressed by its standard deviation, a, and the 
number of years of observation, N, on the assumption that the observations 
are random events following the normal Gaussian frequency distribution. This 
course was adopted because many natural phenomena have frequency dis- 
tributions which are approximately normal if the order of occurrence is dis- 
regarded, and it was thought that the examination of this case, which is a first 
approximation to those occurring in nature, would indicate the form of relation 
to be sought. Gauss’ normal frequency curve (which is also the curve for the 

| frequency distribution of errors of observation) gives the frequency of depar- 
tures of a quantity from its arithmetic mean in relation to their magnitudes. 
If y is the number of departures lying between x + 3 andz — 3 when the total 
| number of values of the quantity is N, the equation of the curve is 


ye (No N 2m) 6-8. Scat see pam (1) 


Incidentally, the assumption that the discharge of a river is a random event 
is the basis of the works of Messrs. Hazen‘ and Sudler;* under this same assump- 
tion a theoretical solution has been found. 


3. Tus Frequency DistripuTions oF RANDOM EVENTS 


It is well known that the Gaussian normal frequency curve (Eq. 1) fits 
' many distributions of events in which success or failure are equally likely and 
in which each event is independent of the next—for example, the frequency of 
accidental errors in a set of observations. It is also well known that an approxi- 
mation to this curve is given by the terms of a binomial expansion. When a 
set of m coins is tossed, the probability of the occurrence of r heads and (m — r)- 
tails at any throw is mC; (3)’ (2)""” in which »,C, is the number of combinations 
of m different things taken r at a time. 
If the set is tossed N times (the number NV being fairly large), the average 
frequency of occurrence of 0 heads and m tails, 1 head and (m — 1)-tails, etc., 
is given by the terms of the binomial expansion: 


NALD =N A+ mCi t wC2 + 0+ + Cm + 1) CH)" vse (2) 


8 RESERVOIR CAPACITY 


As m increases, Eq. 2 approximates the normal curve, and if m = 10 the ap- 
proximation is fairly good; in fact, for many purposes it is good enough. 

Table 3 shows how the number of heads would be expected to be distributed 
on the average when a set of 10 coins is tossed 1,024 times—that is, the binomial 
distribution. The distribution derived from the normal curve, and one actually 
obtained by trial are also given. 


TABLE 3.—CoMPARISON OF FREQUENCIES 
SS —————————————————————————————_——_—_—__ TEE 
NuMBER oF HEaps, r 
Distribution 


Binomial. ...........-- 1 10 45 | 120 | 210 | 252 |} 210 | 120 | 45 10 1 
Normal Gaussian....... 2 10 43 116 | 212 | 258 | 212 | 116 | 43 10 2 1,024 
Actual experiment...... 0 16 53 | 122\| 209 | 240 | 194 | 138 | 46 5 1 


In the normal distribution the numerals represent the nearest integer. The 
relation between the binomial and the normal distribution has been treated 
by G. Udny Yule,’ who defines the normal distribution as the limit of the 
binomial distribution when m is indefinitely increased. For theoretical pur- 
poses the computer can use whichever form is the more convenient. 


4. CALCULATION OF THE RANGE OF CONTINUED SUMS OF DEPARTURES 
oR STORAGE FOR A NoRMAL FREQUENCY DISTRIBUTION 


Assume a number of consecutive values of an element whose variation from 
its mean is distributed normally, as in the binomial expansion already dis- 
cussed, but where there is no correlation between successive values or groups 
of values, as in tossing coins. For example, consider the definite case of tossing 
a set of 2 m coins N times, calling a head a gain and a tail a loss. For present 
purposes each toss of the 2 m coins is recorded as the number of heads minus 
the number of tails, and the departures of this quantity from its mean are 
added to find the range R between the maximum and the minimum of this 
summation curve. The range thus found, from WN tosses of a set of 2 m coins, 
may differ slightly from the range found by the summation of heads minus 
tails from 2 N m tosses of a single coin. This can be easily seen by tossing a 
set of 10 coins one after the other, and recording individual results. Then the 
sums of heads minus tails are the same at the end of every set of ten whether 
individual tosses are considered or not. It is as if the continued sum were 
plotted using single tosses, every tenth point being marked. The curve ob- 
tained by taking the tenth points would be a close approximation to the curve 
of individual points, but the latter would usually have a slightly greater range. 
For example, suppose that the maximum of the curve was produced by aset of — 
6 heads and 4 tails which would put the end-point of the set two units above. 
the end-point of the preceding set. If the 6 heads occurred and then the 4 
tails, there would be a point on the curve recording individual sums of heads — 
minus tails which would be six above the end-point of the previous set—that is, 


(ety Gi i eR a= 
“An I sAaaanie - 
od Ed., 191 eomieon to the Theory of Statistics,” by G. Udny Yule, Charles Griffin and Co., London, 
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the range would be increased by four. For the present this difference is 
ignored, but it will be considered later. ; 
In the case selected for analysis 2 m coins are tossed N times, heads and 
tails or gains and losses are equal at the end of the trial, and nis written for N m. 
The number of orders or arrangements in which n gains and n losses can 
occur is 


in which 2,C, is the number of combinations of 2 7 different things taken n at 
a time. Among these orders there are 2nCay, instances in which, at some 
point in the process of tossing, losses exceed gains by / or more.’ 

There are three cases: 


_Case Condition No. of arrangements 
: . anln 
a Losses never exceed gains oe 
(a) . n+1 
: ae! 
(b) Gains never exceed losses a 
n+l 
(c) Losses never exceed gains C 
° Qnvan 
and gains never exceed an— 1) 
n+1 
losses 


In case (a) the range is the maximum value of gains minus losses; in case (5) 
the range is the maximum value of losses minus gains; and in case (c) the 
range is the sum of the maximum of gains minus losses, and the maximum of | 
losses minus gains. 

Summing all values of the range to obtain a mean value, and writing G for 
number of gains and L for number of losses, 


maximum (G — L) + maximum (L — G) 


a2nvn 


Mean range = 


a 2 maximum (G — L) (4) 


QnvUn 


To find the mean range see Table 4. There are 2nCn41 arrangements in 
which, at some point, gains exceed losses by one or more; and there are onCn42_ 


TABLE 4.—CompuraTIONsS TO DETERMINE THE Mean RANGE 
ES re TT ET 


Gains minus losses No. of arrangements Products 
1 onCn+1 = onCn+2 1 (onCn41 = anCn+2) 
2 onCn+2 —~ nvn+s 2 (anCn42 =F onCn+3 
3 anCn+3 — nUn+4s 3 (anCn+a = anCn+4) 
mil eC snct — onCin (n — 1) (onCon—1 — onCon) 
n 2nU 2n TN (2nU 2n 
ee enn 
Sum . onC'n41 anCn41 + anCn+2 + Cie + onCon 


ee a 
8 “Choice and Chance,” by W. A. Whitworth, Cambridge, England, 4th Ed., 1886, Chapter on Priority. 
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in which gains exceed losses by two or more. That is, there are onCn41 — onCn42 
arrangements in which gains exceed losses by one. ® 

Referring to Table 4, 22> = 1 + anC1 + onC2 + * ++ + anCn-1 + anCn + onCn41 
4 +. + Con. Hence, the sum of the products is equal to 3 (2! — onCn). 
The arrangements in Table 4 include all except those in which gains never 


. On . . 
exceed losses, whose number is ; +1 (or onCn — 2nCn41), and in which the 


maximum of gains minus losses is 0. Hence, the mean maximum value of 
Q2n-1 tes te & iC 
2 2nvn 


2nUn 


gains minus losses is and this is also the mean maximum of losses 


minus gains. Therefore, the mean range is 


Since n is large, Eq. 5 can be simplified by the approximation for n! of James 
Stirling which is as follows: 


nl = Vine (;) eer eda 2? (6) 
Hence, 
Ors (2n)! ss V4n 7m (2 n/e)” is 22n (7) 
2nvUn n! n! [Vinx (n/e)"? a oe wis) (0) ye hiwarel Se ee 
and 
R=Vn07 —l=VN mace eee (8) 


since the product N m is large. 

Thus the average range of the accumulated sums (number of heads minus 
number of tails), when 2 m coins are tossed N times (N m being large), in- 
creases as the square root of the number of tosses, exactly like the accumulated — 
error on a line of leveling. | 

The standard deviation of the binomial distribution produced by tossing 
2 m coins is 


: 
This is the standard deviation of the number of heads (r) or tails (2 m — r). 

The standard deviation og of the number of heads minus the number of tails is _ 

twice g,,or V2m. Substituting for V2 m in Kq. 8, 


R=oa Viner =1250aVQN....... pose tey (10) 


In Eq. 10, # is the range of the continued sum of 2 m N tosses of a single 
coin, and is a little larger than the range obtained from N tosses of a set of 
2m coins. The average difference d between the two computations for range 
can be found as follows: 

Toss a coin a large number of times and plot the continued sums of heads 
minus tails, marking a given number of them off in sets of 2m. The set of _ 
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2m that contains the maximum sum is selected for scrutiny. If there are r 
heads and t tails the end-point of the set is r — ¢ higher than the end-point of 
the ee set. The number of orders in which r heads and ¢ tails can occur 
IS amUr. 

The number of orders in which at some stage of the tossing the number of 
heads shall exceed the number of tails by A or more is omCmin. At the end of 
the set, heads exceed tails by r —t = 2r—2m. The number of orders in 
which at some stage heads exceed tails by more than 2r — 2 m is 2mC'2r—2m41- 
In all these cases the maximum of the summation curve of individual tosses 
will exceed the maximum of the curve derived from sets of 2 m tosses. The 
/ mean value of the excess d can be obtained by a procedure similar to that 
already employed in finding the mean range of the summation curve. The 
result of this is 


Peele [eee am) ha (2 m)! 
= ga | Gai lt 3 Gain FM 
2m)! 
$5 Ot + @m—3)2m+Qm—11]..c0) 


There is a similar effect on the minimum of the summation curve, so that 
the average range of the summation curve of individual tosses is greater than 
the range of the curve from sets of tosses by 2d. Some values of 2d are given 
in Table 5. The correction 2d to ; 
the range is independent of N and TABLE 5.—TypicaLt VALUES oF 2d 


so has a less proportionate effect as 


N increases. VALUES OF ™ 
Magnitude 
5. EXPERIMENTS WITH RANDOM ' 4 5 6 10 
EVENTS Ae Rena rea is | uz | 20 | 27 
Correction*..... 4.1 4.3 4.5 4.8 


Table 6(a) gives the results of an 
actual experiment in which sixpences « Percentage correction to R when N = 100. 
were shaken in a box and thrown = ———————————____—_—_—_—_———— 
onto a table. They have been ana- 
lyzed in some detail to give a practical example of the foregoing theory. 
Their frequency distribution has been given in Table 3, in which there were 
1,024 throws of 10 coins (m = 5) which resulted in a distribution close to the 
theoretical. 

The observations were recorded as they occurred and the continued sums 
were computed from 0 to 1,000 observations. In this case (because, at the 
finish, heads and tails are nearly equal) the actual values (heads minus tails) 
are taken as departures, and not the departures from the mean (Col. 1, Table 
6(a)). At the finish there are 74 more heads than tails; so the mean value of 
h — t per toss is 0.074. This is corrected graphically by the aforementioned 
method. Fig. 6, discussed subsequently, shows the result of one of the 
experiments. 

In each case R (Cols. 3, Table 6) is corrected to what it would have been if 
the departure had been reckoned from the mean. The theoretical value of o 
is 3.16 for a binomial distribution. 


i 
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Table 6(a) shows that R/(o VN) tends to be constant and the values from 
1,000 trials (Cols. 4) agree closely. with the theoretical values (Cols. 5). The 
trials illustrate a point not examined in the theoretical work, and that is the 
variability of o, R, and R/(c VN) from set to set of 100 throws. 

The maximum variations, as percentage departures from the mean, are as 
follows: o, 11%; R, 50%; and R/(o VN), 39%. The standard deviation (or 
standard error) of ¢ for a set of 100 throws is about 7%. 

Another device for obtaining a normal frequency distribution is what may 
be called a probability pack of cards (Table 6(b)). In this pack the cards are 
numbered +1, —1, +3, —3, +5, —5, +7, —7, +9, and —9, and the numbers 


TABLE 6.—EXPERIMENTS WITH RaNpoM EVENTS 


DEFINITIONS OF COLUMN HEADINGS 
Cols. 1. Primary Observations: 


Table 6(a), Excess of heads minus tails = h —t. 
Table 6(b), Sum of recorded numbers = 2. 
Table 6(c), Excess of even numbers over odd numbers = H. 


Col. 2. Standard deviation of corresponding primary observations (Col. 1) =o. 
Col. 3. Range of accumulated deviations = FR. 


Col. 4. Computation of wid 


from observed data. 


af 4 
Col. 5. Theoretical values of data in Cols. 4 corrected for the effect of tossing in sets. 


(a) Ten Stxpences Tossep | (b) Prospasiuity Carps Cur (c) Szrra, NuMBERS 
1,000 Times 1,000 Times oF Bonps 
No. of Dees eg GE Bae Se OT RET DET RRTU EY Fre ETE 
Trials R N R N R/(oVN 
hod e AR ce Za) el) BOL ee | ae oe 
(1) | (2) |. (3) | (4) |) | @) | (2) | (3) | (4) | &) | G@) | (2) | 3) | @ 4] ©) 
1-100 26] 3.08] 39 1,28] 1.20] —26] 3.16] 40 1.27| 1.20] —10| 3.42) 39 1.14} 1.20 
101-200 24{ 3.10} 25 0.81 -| +38] 3.00} 30 1.00} . 52 | 3.22] 32 0.99 
201-300 —26| 3.36) 44 1.31 - 3.99] 35 0.88} . — 4} 2.30) 24 1.04 
301-400 —30 | 2.88) 29 1.01 —44| 4.06] 36 0.89} . —30 | 2.86) 50 1.75 4, 
401-500 +24 | 3.59) 40 5 A .| +28] 3.87] 65 1.68] . 18 | 2.79| 47 1.68 
501-600 —32 | 3.06} 31 1.01 _ 3.58] 48 1.34] . 6} 3.13] 56 1.79 
601-700 +12] 3.56) 44.5] 1.25 .| +10] 3.57] 33 0.92} . 8 | 2.97) 46 1.55 
701-800 +42} 3.02] 28.5] 0.94 —22| 3.11] 44 1.42] . 30 | 2.78) 44 1.58 
801-900 38 | 3.03] 25.5 | 0.84 .| + 4] 3.82] 78 2.04] . —48} 3.02) 49 1.62 
901-1,000 | — 4] 3.50} 54 1.54 —14| 4.25) 48 1.13] . 8| 3.44! 40 1.16 
Mean 2. | 3.22] 36.1 1.11] 1.20] ....] 3.64] 45.7] 1.26} 1.20 2.99} 42.7) 1.43] 1.20 
1-500 +18 | 3.20} 87 1.21} 1.23] — 6] 3.62] 81 1.00} 1.23 26 | 2.95] 84 1.27} 1.23 
501-1,000 | +56 | 3.23] 94 1.30] 1.23] —24] 3.67] 111 1.35] 1.23 4| 3.09] 87 1.26} 1.23 
1-1,000 | +74] 3.22) 121 1.19] 1.23} —30} 3.64] 116 1.01} 1.23 30} 3.02] 88 0.92] 1.23 | 
Mean obeys Ve aieae Tee tote yf MeL Oy DBO ey ce Praag, Wek eesty ed}, LAS re esol fg eacd eet tee OND aoa ‘ 


_ of each kind are proportional to the corresponding ordinates of a normal fre-. 
quency curve. They are: 13 one’s; 10 three’s; and 6, 3, and 1, respectively, 
of the others. The approximation of these numbers is fairly close. The cards 
are first well shuffled and then cut, and the number on the exposed card is 
recorded. The cards are reshuffled slightly and cut again, and so on. The 
numbers recorded may be taken as corresponding to observations of a quantity 
whose frequency distribution conforms to the normal Gaussian curve. 

This is a quicker process than tossing and counting coins. To toss 10 coins 
100 times required about 35 min, whereas shuffling and cutting for 100 cards 
required 20 min. Table 6(b) records the results of 1 trial of cutting 1,000 times. 
As in Col. 3, Table 6(a), in each case R is corrected to what it would have been 


é 
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if the departures had been reckoned from the mean; and, again, R/(o ¥ N) in 
Cols. 4 approximates the value obtained theoretically for the case when N is 


large (Col. 5). 


A third set of chance events was examined (Table 6(c)). From time to 
time newspapers announce the numbers of bonds that have been drawn for 
redemption. A collection of unconnected random events was made by dividing 
the numbers of the bonds into groups of ten, as follows: 


(b) 


fon 
“ 
S 


PP Ph Pp PP RB CO CO OO 
1 ONOORNHF OOME 
LwNOntoawornooworn 
LNworrhanooos 
LPP PP PP & WK! OO 
eR a. ee 
NONWHOAHONOS 
1LeOHADOWH DW 


Number of Even Digits Minus Number of Odd Digits 
—2 +0 +2 2 -—-4 —2 


The left-hand column of a group was ignored as only a few different digits 
are likely to occur in it. The number of even digits less the number of odd 


‘digits in columns (a) (0) (c), (d) (e) (f), etce., are recorded. This difference 


yields a sequence of numbers exactly analogous to the sequence previously 


- obtained by tossing coins, and the labor of producing it is much less. 


By this method a sequence of 1,000 random numbers composed of the even 
numbers from 0 to -+10 was made and this was used as before. Table 6(c) | 
gives the results of this trial. 

In each case R (Col. 3) has been corrected to what it would have been if 


_ the departures had been reckoned from the mean. The theoretical value of 


R/(c VN) (Col. 5) has been corrected for the effect of tossing in sets. 
Although the mean value of R/(¢ 4 N) found from the three sets of random 
events in Table 6 is very close to the value found by mathematical analysis 
(Cols. 5 compared to Cols. 4) yet individual values for sets of 100 observations 
vary considerably. The extreme ranges in the three cases are as follows: 


R 
Table ‘ oVN 
Ga eee i. oa ee coe eee WR eke Shrek 0.81 to 1.31 
GUC) ay eck avis s PAtapeak hom Dike he we ae 0.88 to 2.04 
(7 SULA Panera eta a ag Pan 0.99 to 1.79 


The mean value of R/(¢ + N) from the thirty sets of 100 observations is 1.27 
and the standard deviation is 0.32. The relation between R, o, and N for 


random events having been established theoretically and having been found 
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to agree with practice, it is now necessary to find in what respects natural 
phenomena are similar to random events and in, what respects they differ. 


6. CHARACTERISTICS OF NATURAL PHENOMENA 


The phenomena that chiefly concern the subject of this paper are river 
discharges. Unfortunately records of these rarely cover periods of more than 
50 years, although levels (which are usually an indication of the discharge) have 
been recorded over longer periods. The discharge of a river is closely related 
to the rainfall on its basin, and rainfall records for many places cover more than 
100 years. The investigation began with these phenomena. 

The main similarity of river discharges, levels and rainfall, and random 
events is that their frequency distributions conform approximately to the 
normal curve, although they are usually slightly skew. As a rule the observa- 
tions are too few to give very regular 
distributions, but the flood levels 
(nearly 1,000 of them), of the Nile 
River at Cairo, Egypt, recorded on 
the Roda gage, give a very regular 
distribution which is well fitted by 
a normal curve. This distribution 
is shown in Fig. 2. 

Although in random events 
groups of high or low values do 

So eT SS vcour, their tendency to occur in 
Difference from mean in metres, natural events is greater. This is 

Fie. 2.—Maximum ANNUAL Gace Reapina ar _— the main difference between natural . 
SE gn ait SION ee and random events, an example of 

1946 A.D. which is the discharge of the Nile 

River at Aswan. The long series of 

records of flood levels at Cairo shows the same phenomenon.® There is no 
obvious periodicity, but there are long stretches when the floods are generally 
high, and others when they are generally low. These stretches occur without 

any regularity either in their time of occurrence or duration. 

In the investigations of Messrs. Hazen and Sudler the frequency character- | 
istics of river discharges are assumed to be like those of random events, and — 
this is the common assumption in estimating probability of high or low floods. | 
It will be shown that this assumption must be used with caution and, in cases in 
which storage over long periods is concerned, it is only an approximation. 


Frequency per thousand. 


7. THe Rance (R) oF THE ContTINvED Sums oF DEPARTURES 
FoR NaturaAL PHENOMENA 


The investigation began experimentally, before any attempts to find a 
mathematical theory were made, by finding R for any available long series of 
river discharges. Such series were scarce and so the work was extended to. 
rainfall data, which are more plentiful. When it was found that rainfall data 
en ag eS ee 


9 ee Records of the River Nile,” 


p. 1012. by C. 8. Jarvis, Transactions, ASCE, Vol. 101, 1936, 
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- gave results similar to the river data the work was extended to other long series 
of natural events. 

For example, Fig. 3 shows a summation curve of departures for the maxi- 
mum levels of a lake on the Dalalven River in-Sweden. There are 176 years 
of observations, divided into periods of 44 years, 88 years, and 132 years; and 
the ranges, R, of the summation curves for the departures from the mean in 
each of these periods are found. 

The sums of departures from the base OX are plotted.. For the sum of, 
departures from the mean for the entire period the curve is referred to the 


_ 1000 


500 


»¢ 


O By 
1764 7800 1850 7900 1940 
(meee ey see Sat! ee ee ee 
oO 44 &8 132 176 Years 


Fig. 3.—Summation Curve or Departures, Maximum LEvELs OF A SMALL 
LAKE ON THE DALALVEN River IN SWEDEN 


| axisOD. Similarly, the summation for the first 44 years is referred to line OA, 
for the second 44 years to line AB, for the first 88 years to line OB, and so on. 
The range, R, of the last period BD is QN + PM. 


8. VARIATION OF R witH N AND o 


The theoretical investigation shows that R is proportional to o in the case 
of random events, and this must be approximately the case for any events that 
have normal or nearly normal frequency distributions. The device of express- 

‘ing deviations from the mean in terms of o is a common one, derived from the 
“expression for the normal curve; and, since R is a function of these deviations, 
R/o has been calculated. Complete tabular records of the computation in- 


cluding mean values of N and K (- eee) are in the complete original 
manuscript, cataloged for reference in the Engineering Societies Library.” 
Seventy-five different phenomena were used, and 690 values of R were com- 
: puted for different values of NV. 


: 
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A summary of the results is given in Table 7 in which the phenomena are 
divided into groups according to their nature and the value of N, and log (R/c) 
is plotted against log N in Fig. 4. 

In general, R has not been computed for periods of less than 30 years be-. 
cause, for very short periods, it has little practical value although it may have 
some formal value. For N = 2, R/o has the exact value, unity. The longest 
period used is 2,000 years and this is for clay layers in Lake Saki in the Crimea. 
Fig. 4, in which the means of small groups of observations of the same kind 
and length of period have been plotted, shows that log (R/c) is a linear func- 


TABLE 7.—Rewation oF R/c aNp N For GRouPs OF PHENOMENA 


Phenomena Newt eal — Log N Log — K 
(a) GrouPp or 99 Caszes 
«1 8] S| | as | 8B | oe 
i Is, discha , and runoff....... 8 K * i i 
ia rite tae 8 62 13.1 1.79 1.08 0.72 
: 108 16.4 2.02 1.19 0.69 
: 12 105 19.6 2.02 1.27 0.75 
10 208 36.5 2.32 1.54 0.77 
9 309 53.9 2.49 1.72 0.79 
8 420 60.3 2.56 1.77 0.78 
Roda nilometer.....5.......2cc cece eee (é 511 67.7 2.71 1.82 0.75 
6 613 81.3 2.79 1.89 0.77 
5 716 104 2.85 2.01 0.79 
4 820 122 2.91 ' 2.08 0.80 
3 927 129 2.96 2.11 0.79 
2 1,040 130 3.02 2.12 0.78 
Mean of 99 cases.......... 02sec ee eee vests J onoke eee 2.28 1.48 0.75 
(b) Group or 168 Caszs 
Rainfall stations with one value of R-. ... . 7 71 14,0 1.82 1.12 0.74 
Two groups of values of R®....... sake Hee i = ye ae he Ki pe 
Rainfall stations with three groups of | (66 38 8.2 1.57 0.91 0.72 
Waliies Ol seine ce tetaie ieicle Soak alse 38 78 14.4 1.87 1.13 0.72 
25 121 22.2 2.08 1.31 0.74 
Mean of 168 cases................2.-.04. cra are 1.79 1.08 0.70 
Sade cheat is BRAK SOE DT ghee SARS BRS ee eebet Be Se A GEREN aR a SS eT 
(c) Group or 109 Caszs 
Ee oe a Oe 8 Ee a Se i ee 
Temperature and pressure stations with 3 27 37 7A 1.55 0.85 0.68 
VOIUGR OL TE Lass Ca sera ioe ois Se eal ane hee 18 73 12.0 1.86 1,06 0.68 
at ” ice es nae 0.71 
Temperature stations with 4 values of R. . 12 88 15.1 1.94 V7 Ot 
uy i aoe oh 1.32 0.72 
Mean of 109 cases.............0002000ee rie aes 1.81 108 070 
ANN 


# | 8] Be] a8 | re | Be 
Annual growth of tree rings............. 14 200 45.0 2.30 Ve O82 
4 300 79.9 2.48 1.90 0.87 
4 462 118.6 2.66 1.87 0.79 
Mean of 85 cases................000c005 ere as Rises Re io 736 O00 
Ht 


“¥ 
Fy 
EY, 


> 


RESERVOIR CAPACITY 17 


TABLE 7.—(Continued) 


No. of N 
Phenomena akon years = Log N Log— | K 
(e) Group or 90 Cass 
; ; 44 50 10.9 1.70 1.02 0.73 
Thickness of annual layers of mud; Tamis- 22 100 21.3 2.00 1.31 0.77 
kaming, Ont., Canada, and Moen, in the 11 200 42.7 2.30 1.58 0.79 
Sogne District, Norway.............-. 4 300 82.5 2.48 1.90 0.87 
3 400 126 2.60 2.09 0.91 
4 550 115 2.74 2.00 0.82 
4 1,100 181 3.04 2.19 0.80 
MWrean of 90 casess.c. 6. theese cele sae es si 1.98 1.30 0.77 
sake SN Me Aa ce eM Ws NP RS PR SOM, PRN Se tS AE Eat 
(f) Group or 114 Caszs 
Eat Soles RCI EIN Nabe Mc TR SE i Te es snes ee sk 
40 50 9.7 1.70 0.98 0.70 
Thickness of annual layers of mud, Lake 40 100 - 15.3 2.00 TZ: 0.69 
Saki in the Crimea..............-...- 20 200 25.0 2.30 1.39 0.70 
8 500 47.9 2.70 1.66 0.69 
4 1,000 84.0 3.00 1.91 0.71 
2 2,000 179.0 3.30 2.24 0.75 
MMfean- of 114 casesiii fb icie oe tle tees es Sahel Savas Se 2.06 1.22 0.69 
SS Ee a i, on ret re [eee ee eS ns 
(g) Group oF 25 Cases 
ems een ee ee ee 
12 64 1 ay iep 1.06 
~ Sunspot numbers and wheat prices....... 6 124 22.1 2.09 1.34 0.75 
te 237 16.9 2.36 1.22 0.60 
Mean of 25 cases.........----0+:- cece eee sar 2.0 
See — — — —— = nn I cn = lt Ee 
(h) Group oF 259 Casxs 
ete ee rs a a eee 
Weight mean of 259 cases*.......-.--.-. | 97 17.8 | Be a | 


Cee eee eee 


« Rainfall stations with one value of R; includes temperature at one station. © Rainfall stations with 
two groups of values of R; includes temperature and one pressure. ¢ N ranges from 81 to 120. m 


r ee ——————_______ 


tion of log N and can be expressed as 
R 
log me K log 


In each main group of observations a line drawn through the point log (R/a) =0, 
log N = 0.30 (N = 2), and the center of gravity of all the observations is a 
good fit. The slopes of these lines (K) vary from 0.69 to 0.80. 

Whether there is any theoretical significance in the fact that K is approxi- 
mately 3 is not known. The value of K has been calculated for each of the 
690 individual values of R/o and the frequency curve shown in Fig. 5 has been 
drawn. The data for this curve are as follows: 


Description K 
Number of values...........--2 eee ee eeeees 690 
MV GARAIIIO | th tie Os Rk tees ieee.) ks a) Vie mes 0.729 
Standard deviation.............--2+eeeeeeee 0.092 


Range of individual values...........---++--: 0.46 to 0.96 
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Lo3(R/ao) 


/ 2 3 log N 


2 5 105720 50 100 200 500 1000 2000 


Fig. 4.—ReELation BETwEEn RANGE oF SuMMATION CURVE R, THE STANDARD DEVIATION o, 
AND THE YEARS OF RecorpD N 


Fig. 5 shows that there is a slight skewness in the distribution of K, the 
mode tending to be slightly greater than the mean. However, it will be seen 
that the normal curve which has been fitted is a close approximation to the 


observed distribution. Referring to Table 7, a summary of K-values can be 
compiled as in Table 8. 
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Table 8 and Fig. 5 show a remarkable result—that the values of K are 
nearly the same regardless of the class of phenomenon from which they are 
determined. It is necessary to examine them a little further to see if the differ- 
ences are significant. It will be noticed that the standard deviation of K from 
all the values is a little larger than 
that for most of the values for the 100 
groups. When the standard devia- 


tion of the means of groups is cal- 8 
culated it is found to be about 4.25 8 
times the standard deviation of the S 
mean of a similar group of random & 
values of K. This suggests that K @ 50 
does vary slightly with different 3 
phenomena. The largest mean value z 


of K is 0.81 from tree rings, and 
this is derived from four sets of trees, 
each set being from one of four dif- 
ferent places. The means of sets : : a 
are 0.84, 0.80, 0.81, and 0.81, so K 
that it would appear that the de- Fic. 5.—Frequancy or Inpux K 
parture of the trees from the mean 
of all the phenomena is significant. Similarly, the difference in K between 
the varves of Lake Saki (the Crimea) and those of Tamiskaming (Ontario, 
Canada) and Moen (Norway) is greater than would be expected on a chance 
distribution. Tamiskaming and Moen are significantly different. The rainfall 


TABLE 8.—Summary oF K-VALUES 


Standard 

Phenomena ; No. Mean aevintionl 
River levels, discharges, et... ...---. esse reer es ree erste resets 99 0.75 0.077 
TRY HSL voice Oe BIB DECI bid SOM tne cick ots -aU Gana Meir aaa 168 0.70 0.069 
Temperature and pressure......---++eseeeeseerrtrsestrttts Betas 115 0.70 0.085 
Annual growth of tree rings....-.---+-+eee reer er cers este 85 0.81 0.078 
Varves (Lake Saki in the Crimea). } ob. ccc le ns nee ete cee e eee dns 114 0.69 0.064 
Varves (Tamiskaming, Ont., Canada, and Moen, Norway)...-.--++ 90 0.77 0.094 

Sunspot numbers and wheat prices (combined as miscellaneous 

Plier onic a) eave iere ors aisles sisal ete = cinta ieee Sb peal ts 20s 25 0.69 0.086 
PNReR Stand tOtalsee a als cise eywiti=i¢ vies 20s ccerasnicneoderereretay= sien “iaisie/ ete) 4p, 210 690 (0.729) (0.092) 


ee ao 


observations have been examined in detail, with the result that the means for 
fifteen stations with some 7 values of K to each station do not differ greatly 
from a random distribution. 
The general conclusion is that, if a long series of several hundred years is 
available for a phenomenon, it is probably best to take the mean value of K 
from these as an indication for the future. However, if only a short period of 
from 100 years to 200 years or less is available, the best plan is to take the mean 
value for K from all available material. ‘The mean adopted for this purpose 
is 0.72 instead of 0.73, in order to give a little more weight to the more precise 
measurements relating to rivers, rainfall, and temperature. The mean value 
adopted, therefore, is K = 0.72 + 0.006. In using this result it is to be noted 


vy ae 
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that it is the mean value of a quantity that has ranged from 0.46 to 0.96, and, 
therefore, some margin of safety must be allowed. om 

The probable error of a single value of R/o can be derived by differentiating 
Eq. 12 with respect to K and using the data in Table 8, which yield a probable 
error of 24%. Confirmation of this value is obtained by taking the values of 
R/c for rainfall, river discharge, temperature, and pressure for values of N 
between 80 years and 120 years, and correcting them to a period of 100 years. 
The probable error calculated from these corrected values is 28%. Fig. 10, 
introduced subsequently, shows how R/o can be reduced by 23%, by reducing 
the draft by 0.04 c. 

In Fig. 5. the mean representing all the lines is represented by 


2 
which means that R/o increases more rapidly with N in the case of natural 


phenomena than in the case of chance events for which ‘the power of N = 3. 
The value of R/o for N = 100 from Eq. 131s 16.8. In the computations previ- 


ously cited? the equation R/o = 1.65 4 N was used, which gives R/o = 16.5 


for N = 100. 
If a phenomenon is exactly periodic, with constant amplitude, R will have 
a definite maximum which will be reached at the end of each complete period. 


In this case R will not vary in the manner already found for natural phenomena. — 


The phenomena considered in this paper, however, are not of this type. In 
the case of sunspot numbers (which are the nearest approach) the periodic part 
does not have a constant amplitude, and has irregular variations superposed. 
These irregularities cause R to increase as N increases and the mean value of 


0.72 
hin (*) = OGL Nit es ee nae (13) 


K is 0.75, with no indication that the phenomenon is different from the many ° 


others considered. 

The very striking result from the investigation is the small variation of the 
index K over such a wide range of phenomena, and hence the possibility of its 
utility in many fields not mentioned in this paper. 


9. THE VARIATION OF THE MEAN AND THE STANDARD DEVIATIONS 
DERIVED FROM SHORT-TERM RECORDS 


When records of a natural phenomenon extend over long periods there are 
considerable variations of both means and standard deviations from one period 
to another. The primary practical interest of this paper is in water supply 
and, hence, in rainfall. Unfortunately 200 years is about the longest series of 
observations obtainable for rainfall, and few stations reach 150. years. Some 
of these records are given in Table 9. They are divided into periods of 40 
years or 50 years to illustrate the recorded variations. This division is chosen 
because quite frequently water supply or hydroelectric projects have to be 
based on such short records since nothing more is available. The reservoirs in 
the Great African Lakes are a case in point. 

The average range of variation of the mean values of the rainfall in Table 9 
- is about 14%, whereas the standard deviations vary on the average by 23%. 


& 


\ 
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‘This is more than would be expected from the ordinary theory for random 
distributions. The same is true in the case of temperature. 

The uniformity in the values of R shown in the preceding section encourages 
the examination of some of the long-period records, a summary of which is 


TABLE 9.—VARIATION OF MEANS AND STANDARD DEVIATIONS 
IN THE Case OF RAINFALL 


Duvetion b- Langth wg MBEAN STANDARD DEVIATION 
Station of records | of period oe a 
(years) (years) Berioes 5 ie : hig 
Maximum | Minimum | Maximum | Minimum 
Philadelphia, Pa........ 126 42 3 42.7 41.4 7.0 5.6 
Boston, Mass........... 128 43 3 44.2 38.2 TA 5.0 
~ Madras, India.......... 133 44 3 51.3 48.4 15.6 14.3 
Rome, Italy.........-.. 151 50 3 886 777 17.1 15.7 
Stockholm, Sweden..... 161 40 4 57.8 38.3 13.5 7.8 
Padua, Ttaly sede aries uel 43 4 940 793 179 147 
Milan; Italy... 25% 0.3). + 171 49 6 104 95 21.2 15.1 
Zwanenbourg, the 
Netherlands.......... 211 53 4 77 69 13.1 10.8 


ede ee a ee 


compiled in Table 10. In Table 10(a) the variability of means for 50 years 
and 100 years is compared with what would be expected if the frequency were 
‘random, and not affected by any tendency for high or low values to be grouped 
together. 
The computed value of the standard deviation (c) of a mean for 50 years 
(Col. 4) is obtained by taking the standard deviation obtained from all the 


TABLE 10.—VarIATIONS OF MEANS AND STANDARD DEVIATIONS 


(a) SranpARD DEVIATIONS (b) SranparD DryIATIONS OF 
or MEans STanpDARD DEVIATIONS 
No. of 50-YEAR 100-Y EAR 50-YnaR 100-YEAR 
Phenomenon years, Means Means PERIODS Prriops 
N : Sree ee Site § Eee ee 
Ac-| Com-| Ra-| Ac-|Com-| Ra- | Ac- |Com-| Ra-| Ac- | Com- Ra- 
tual|puted| tio | tual]/puted] tio | tual puted] tio | tual |puted| tio 
(1) (2) 3)! 4) 16)1@)| © 16} @) | © 16] @ | @ (5) 
Varves: ; 
1000 (F810 |a7 (io |oe | 1s 
rb a 1,00 : { ozalidh ‘Ble 
Lake Saki, Crimea... 1000 92/13 |1.71/16]09 | 1.7 26 10.6 |42)3.1 |0.5 | 65 
950 2.2 | 1.2 1.8 | 2.0 | 0.9 2.3 
Moen, Norway.....- 1,000 4111.9 121135114 | 25/56 {12 | 48/45 0.9. | 5.0 
Tamiskaming, Ont., 
Canada. is alee 1,200 4210.8 | 5.5 13.5 6 | 6.2 |1.11 |0.36 | 3.1/1.6 |0.3 | 5.4 
Ring Thickness: 
Beanoia trees...... 900 0.12] 0.04 | 2.8 | 0.10] 0.02 | 5.6 | 0.034] 0.014} 2.4 | 0.032} 0.011 2.9 
Binpeiadl pines... :<. 500 0:09) 0104.) 220 fon ed ee si ...| 0.040} 0.029] 1.4] ....]---- |--+- 
Roda gage (free o 
f trend) 3 x pe re es 1,050 0.24] 0.10 | 2.4 |0.21| 0.07 | 3.0 ]....]....|....| 0.145] 0.042 3.4 
INDOBTS Bie rel cicis 2 ee es ieinte Pe All) op 1 (OF eon ol Rett Baa GEO lan eae eS GRP A tPA incre (4.6) 
SE ee ee eee 


individual observations for the period given in Col. 3 and dividing it by V 50. 
Cols. 5, Table 10, are the ratios Cols. 3/Cols. 4. It appears from Table 10(a) 
that the means are much more variable than they would be in the case of a 
random distribution, the standard deviation of a 50-year mean being 2.5 times 
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as great as it would be for a random distribution. Similarly, the standard 
deviation of a 100-year mean is 3.2 times as great as in a random distribution. 
This implies that, over a long period, the maximum value is likely to be higher 
and the minimum lower than would be predicted from the application of the 
ordinary theory of probability to the records from a short period. 


Table 10(b) shows the variability of the standard deviation for periods of 


50 years and 100 years in those cases where the records cover a long time. If 
a distribution is normal and randomly distributed as in chance events, the 
standard deviation of the standard deviation is o/¥2N. This is the value 
called “computed” (Cols. 4). The actual value (Cols. 3) is found in the usual 
way from the standard deviations for all the 50-year periods, etc., and Cols. 5 
contain the ratios, Cols. 3/Cols. 4. Again it appears that the standard devia- 


tion for these natural phenomena is more variable than would be the case in 


a random distribution. 

Summary.—Although many natural phenomena have a nearly normal fre- 
quency distribution, this is only the case when their order of occurrence is not 
considered. Consequently, their representation by a random distribution in 
the usual way is only a first approximation. The tendency to occur in groups 
makes both the mean and the standard deviation computed from short periods 
of years more variable than in the case of random distributions. Thus the odds 
against the occurrence of some particular extreme value (say, a high flood) are 
likely to be less than those calculated from a short period of records on normal 
probability theory. 

It was thought that it might be possible to derive an expression for the 
effect of grouping, involving the correlation coefficient between successive 
values of the variable. Experiments with the natural number 1, 2, ---, N and 

_—1, —2,---, —N, however, showed that this was not the case. The greatest 
value of R is obtained by placing the positive numbers in one group and the 
negative in another and is independent of the order in the groups. The correla- 
tion between successive numbers depends on the order within the groups. In 
the foregoing example, the coefficient of correlation may vary from 1 to 3 while 
& remains constant at N(N + 1). The numbers can also be arranged to give 
small values of R with correlation coefficients varying from 3 to —1. 


10. SroraGE TO GUARANTEE A Drarr Less THAN 
THE MEAN DISCHARGE 


If the discharge from a reservoir over a period is kept steady at the mean, 
the range of the curve of progressive sums of departures from the mean is taken 


as the storage required to keep this constant discharge. The final sum of the 


departures is zero. If departures are taken from a base less than the mean, this 
corresponds to the case in which the reservoir discharge is not allowed to fall 
below this base value. 

The storage required in this case is the greatest accumulated deficit. This 
is illustrated in Fig. 1, of Lake Albert discharges, in which the summation 
curve, or curve of accumulated departures from base 24, is plotted, axis 00, 


Fig. x sae the summation curve for departures from the mean is referred to 
axis OX. 


> ey, tay do 
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If departures are computed from a new base 0.2¢ below the mean, as 
explained previously, the summation curve must be referred to the inclined 
axis OX, where the ordinate XX, is 43 X 0.2 ¢ or 63 milliards. It is obvious 
that, as the inclination of the axis OX; changes, different points on the summa- 
tion curve may become maxima, or minima, and there is no simple connection 
between the storages required to guarantee the mean and to guarantee some 
lesser discharge. This is clearly shown in Fig. 6, in which the result of tossing 


80 Number of Throws 


20 


-20 


20 40 60 xe) 


Fia. 6.—SummMation Curve oF Heaps Minus Tarus, In Tosstne 10 Corns 100 Times 


10 coins 100 times is shown in the form of a curve of the continued sums of the 
number of heads minus the number of tails. 

At the end of the 100 tosses, heads exceeded tails by twelve, so that the 
departures of heads minus tails are reckoned from a base 0.12 below the mean. 
This is corrected in Fig. 6 by drawing the axis OX which is the base for the 
summation curve of departures of heads minus tails from their mean. Simi- 
larly, the axes OX; and OX» can be drawn for the summation curves of de- 
partures from the bases 0.1 o and 0.2 ¢ below the mean. The two points aay 
are the maximum and minimum of the summation curve of heads minus tails 
and also of the summation curve of departures of heads minus tails from its 
mean. ‘The range of the first curve is greater than the range of the second. 
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The maximum accumulated deficit with the mean as base is given by the points 
ab, but when a base 0.1 ¢ is selected below the mean the maximum accumulated 
deficit is given by cei; for a base 0.2 ¢ below the mean, it is de,. The storage 


Fia. 7.—ANALYSIS OF TossED Corns, CarDs, AND DrAwinG oF Bonps To DETERMINE A 
Revation Berween Drarr B anp Maximum Dericir S (Szs Tasux 11) 


TABLE 11.—Vatuzs oF S/R ror Given Ratios (M — B)/¢ 
(Means or 10 Experiments or 100 Triaus Eac#) 


Vautues or (M — B)/o 


Media R 
0.1 0.2 0.3 0.4 0.5 0.6 
COATT San Wenge he Sais 36 0.62 0.47 0.40 0.34 0.31 0.26 
BIOS .S ee cas des. 46 0.71 0.56 0.45 0.36 0.29 0.23 
Botdshisn. ts siigek os 43 0.62 0.48 0.39 0.32 0.27 0.23 
MCANS 25.5005 a Sc lane 42 0.65 0.50 0.41 0.34 0.29 0.24 
OUR seen eta cae —0.19 —0.30 —0.39 —0.47 —0.54 —0.62 


required to maintain a draft below the mean can be computed arithmetically; 
but it is extremely laborious and by far the quicker and more reliable method 
is the graphic one. 


TABLE 12.—Re.ation BrETwEEN Drart B at 


ee 


VALUES 
Magnitude 


0.06 | 0.09 | 0.1 | 0.20 | 0.23 | 0.26 0.29 0.30 0.37 0.¢ 
a Ce ee ee ee ee ee 
(a) Rivers, RAINFALL, AND TEMPERATU 


= 


0.56 0.49 0.46 0.42 0.32 
—0.26 | —0.32 —0.37 |. —0,40 | —0.51 


ae ee el Ge ee aed ee en Oe oe ua 
(b)) TREE Rines anp Varves, Ropa Gs 


SSB Te OLA eine ee LALA Re PS EM GA TSr Te 
YR oe Veshces Se echace 0.64] o48] .... . 
Log (8/R).3... 11; | a. | ewe | SOLIS | PO.B2a fh ec a eee naan Oa ee oi 
SA REA EA Sas ANT Ean PIG 
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11. Tar Revation BeTwEEN DRAFT AND STORAGE 
FoR CHANCE EVENTS 


Let the storage required to maintain a draft B less than the mean (M) be S. 
Values of S (the maximum deficit) have been computed from the experiments 
) with coins, probability cards, and drawing of bonds, as described. They have 
| been calculated for every set of 100 experiments and for values of M — B of 
0.1 c, 0.2c, etc. The ratio S/R has been tabulated and mean values are given 
in Table 11. 
The results are also given in Fig. 7(a) in which log S/R is plotted against 
| (M — B)/c, the relation being 
M-—B 


log 5 = — 0.11 — 0.88 pe ea (14a) 


A little better fit is obtained from the equation— 


| which is shown in Fig. 7(6). Neither formula should be used outside the range 
of values from which it was computed. 
Considerable research was devoted to finding a theoretical relation between 
| S and R, but without success. An extension of the method of Section 5 will 
enable the mean range of the summation curve of departures from a base below 
’ the mean to be found, but it does not result in a simple form. This mean range 
is not what is wanted, however; it is the maximum accumulated deficit in the 
/ entire period, or any part of it that is wanted. The latter is a complicated 
function of gains minus losses, for which no algebraic expression has been found. 
An interesting point in the case of coins, cards, etc., is that, if the mean range 
is called Ri when departures are measured from a base less than the mean, 
then the ratio R/R, is approximately equal to S/R. 


orace on Maximum Dericit S (MEAN VALUES) 


— B)/« 
: Magnitude 
).42 | 0.43 | 0.5 0.52 | 0.53 | 0.55 | 0.6 | 0.63 | 0.66 | 0.68 | 0.7 | 0.88 
winty-EigHt Puenomena; N = 96 YEARS) 
.34 | 0.32 0.24] 0.28] 0.24 0.23 | 0.20] 0.20 0.09 S/R 
9.47 | —0.52 —0.64 | —0.58 | —0.62 —0.67 |—0.73 | —0.76 —0.96 log (S/R) 


svEN PHenomENA; N = 430 YEARS) 
Al ( S/R 
0:88 : | log (S/R) 


j . 
| 


26 RESERVOIR’ CAPACITY 


12. Tue RevaTion BETWEEN Drarr AND STORAGE 
FoR NATURAL PHENOMENA 
The results of a number of computations of the storage (S) required to) 
guarantee a minimum draft (B) less than the mean have been compiled. . 
Table 12 gives the mean values of these detailed computations. Curves of 
10 


00 
S/p 20:97 -0:95V(M-B)/o 


o River discharge 
x Rainfall 
A Temperature 


Log(Y/p)= -0:08-1-00(M-BYo 


o River discharge 
x Rainfall 
A Temperature 


-0-5 0-5) 
~ 
= 
. : 
N Os 
qV® 0 
7) ; 05 (M- B)/e 10 7) 0-5 ViM-BY/o 10 


Fie. 8.—ANA.Lysis oF River DiscHarGe, RAINFALL, AND. TEMPERATURE TO DETERMINE 
A ReuaTion Between Drarr B anv Maximum Dericir S§ 


the exponential and square root forms have been fitted to these mean values 
and the results are shown in Figs. 8 and 9. As far as closeness of fit is con- 
cerned, over the range of observations, there is no significant difference be- 
tween one type of curve and the other. At some future time, it may perhaps 


0-0 1-0 


Log(Sp)= - 0:07-123(M Ba 


0 05 (M-B)/t 10 0 O02 O04 O06 08 10 
ViM-B)/o 
Fig. 9.—Anatysis or Trex Rivas, Varves, AND Ropa Gace To DreteRMINE 

A Revation Berween Drart B anp Maximum Dericir S 


_ be possible to decide that one or the other has some theoretical justification, — 


but this has not so far been possible. The results of the curve fitting may be — 
studied in Table 13. The curves that have been fitted to the means of all the 
observations are shown in Fig. 10. , 
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The equations show very little variation among the three groups of phe- 
nomena. For example, if (M — B)/c is 0.25, S/R from the exponential curve 
for each group is 0.47, 0.47, and 0.42 and from the square root curve for each 
‘group is 0.46, 0.49, and 0.42. These results, again, illustrate a striking 
uniformity. 


TABLE 13.—Comparison oF Curves SHowING RELATION BETWEEN 
Drart B aNp StoraGe on Maximum Dericit S 


Fig. No. Media Log (S/R) = S/R = 
tf Cards, coins, and bonds (Eas. 14)....|— 0.11 — 0.88 (M — B)/o | 0.91 — 0.89 V(M — B)/c 
8 Rivers, rainfall, and temperatures. ...|— 0.08 — 1.00 (M — B)/o | 0.97 — 0.95 V(M — B)/o 
9 Tree rings, varves, and Roda gage. ...| — 0.07. — 1.23 (M — B)/c | 0.91 — 0.98 V(M — B)/e 


0.94 — 0.96 ¥(M — B)/o 


Mean of all observations............ — 0.08 — 1.05 (M — B)/o 


Considering the curves based on the means of all the observations, which 
are given in Fig. 10, it will be seen that in spite of the difference in form, the 
foregoing expressions give practically identical curves over nearly all the range. 

The important result shown by 
all the cases is that a small reduction 
in the guaranteed draft from the 


1-0 


0 Log(Yp)= 0:08 -1-05(M-By/o 


| Maximum value (the mean) makes 
a much greater proportional reduc- 
tion in the storage required to main- 
tain it. For example, a reduction 
of 0.1 o diminishes the storage by 


x Sp 20:94-0.96ViM-bYo} 


OS 


about 35%. This fact provides a 
very practical method of applying a 
factor of safety to compensate for ~ 


the natural variations of RF. 


13. CoMBINATION OF STORAGE 
FROM DIFFERENT SOURCES 


05 (M-BYje 10 


Fie. 10.—Meran or Aut OssERvATIONS ANALYZED 


to DETERMINE A RELATION BetwEeEn DraFr 
anp Maximum Dericit S 


The previous theory has dealt 
with the storage of one variable 
supply, as for example, the water flowing out of Lake Albert. It was shown 
(see Eq. 13) that, on the average, the capacity required to produce a steady 


discharge equal to the mean over a period of N years, on a stream whose 
ae Fe 


"o 
two streams are involved, with variations whose discharges have a normal 
distribution, then if o; and o are their standard deviations, and o is the 
standard deviation of their combined flow, it can be shown that 


standard deviation of annual discharge is g, is given by R =o ( 


a2 = oy + o2g $271,201 02... eee ee ee eee 


in which 1j,2 is the coefficient of correlation between the two discharges. It 
follows from Eq. 15 that, for a given value of N, 


R? => R*; + R*, - 2 T1,2 Ry R, NT Sua Laka sehe: opabauncrohdel otal 
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If the two tributaries vary exactly together, 71. = 1 and R = Ry + Ri. ; 
Except in this case the necessary capacity to equalize the discharge of the main | 
stream is less when it is all in one reservoir on the main stream than when it is} 
in reservoirs on each of the tributaries. It is possible under some circum- - 
stances for the discharge of the main stream to be equalized by a reservoir on | 
a tributary. 

14. PracticaL APPLICATIONS 

As mentioned in Section 1 the investigation arose out of the need to control | 
the flow of the White Nile by reservoirs in Lake Albert and Lake Victoria. | 
The scheme favored by Egypt was to construct a single large reservoir in Lake ; 
Albert of such a capacity as to equalize the outflow over a century. Such a. 
scheme requires the simplest type of regulation, but it was not acceptable to | 
Uganda because of the large area of land which would be flooded for long | 
periods. The Uganda Government proposed the combination of a large storage - 
reservoir in Lake Victoria with a small regulating reservoir in Lake Albert. 
Later a regulator at the exit of Lake Kioga in Uganda was seen to be a necessary 
part of this combination to eliminate the very variable effect. of Lake Kioga on 
the Victoria Nile. The regulation of the combination is more complicated and 
has been the subject of much study. 

In the more simple plan for Lake Albert the capacity of the reservoir would 
be large compared with the natural fluctuations of the lake. The method of 
working would have been to build the reservoir 15 years or 20 years before it 
was required to be in complete action, and, during this time, to fill it partly full, 
leaving a sufficient reserve capacity to control unusual floods that might occur 
after the project was working to full capacity and which would then do damage 
to the Sudan. The simplest working of the reservoir would be to send down a 
constant annual discharge slightly less than the mean, distributing this dis- 
charge as required through the year. This operation would lead to occasional 
years in which a discharge of more than the average would be necessary to 
avoid overfilling the reservoir. 

. This procedure was adopted in the report? for making estimates of the value 
of the projects and of the size of the works required. In these estimates the 
guaranteed discharge was less than the mean discharge by twice the probable 
error of the mean, plus a further small reduction. The actual procedure would 
depend on circumstances, beginning, no doubt, with what was thought to be a _ 
safe draft. If the reservoirs showed signs of reaching high levels, the draft 
would be increased and, if the contents diminished considerably, the draft 
would be reduced. During the early years of regulation, before the full draft 
was required, experience would be gained which would enable a sound scheme 
of regulation to be drawn up based on more information than is now available. _ 

A modification of this simple regulation would include what the writer has 
called “‘virtual storage,” in which an excess on the Ethiopian tributaries can — 
be virtually stored by reducing the outflow from Lake Albert. 

; Lake Victoria, as the main storage reservoir, has the advantage that its area 
is 12.5 times that of Lake Albert. Therefore, a much lower dam will give the 
required storage; but the regulation is complicated by the fact that the natural 
range of rise and fall of the lake is greater than the extra range required for 
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| 
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regulation. This would not matter if there were plenty of room for regulation, 
but the interests of the inhabitants around the lake place a definite limit on 
the height to which the lake can be raised, and the natural range of the lake 


plus the extra effect of regulation must be fitted within the allowable limits. 


This requirement fixes the capacity needed in Lake Albert, where extra storage 
over that required for the annual regulation is necessary to provide against 
unusual floods. The regulation of the two lakes in combination, therefore, 
departs from the simple one of a steady annual discharge out of Lake Victoria 
and a steady annual discharge out of Lake Albert, which would be almost auto- 
matic. The more complicated regulation will depend for success upon a close 
study of rainfall, runoff, discharges of tributaries, and trends of lake levels. 
However, the time that will elapse before the project is working will provide 
additional data, and experience in regulation. 

It is fortunate that the Nile supply comes from two sources—Ethiopia and 


the Great African Lakes—and that, as far as available information goes, the 


correlation between their contributions is small. It is unlikely, therefore, that 
both sources would suffer extended periods of drought together, and this can 
provide a factor of safety against an unusual sequence of low years on the Lake 


Plateau which includes Uganda and parts of Kenya, Tanganyika, and the 


Belgian Congo. To take advantage of this the scheme of projects reported? 
must be considered as a whole and worked as a combination. Should an 
unusual succession of low years on the Lake Plateau cause a serious decline in 
the stored water, it would be necessary to make full use of virtual storage. 
Then, by beginning to store a little earlier in the main Nile reservoirs, the quota 
from the lakes could be reduced and the possibility of a dangerous depletion 
of the storage would be removed. 


15. CoNcLUSION 


This paper has been mainly confined to theoretical matters, and the prac- 
tical outcome therefrom has been reported elsewhere.” 
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APPENDIX. NOTATION 


} 
The following letter symbols have been adopted for use in this paper and in . 
the discussion: 


B = annual draft from storage S on a stream, less than the mean annual 
flow of the stream; 
number of combinations of m different things taken r at a time; 
departure: 
d = average difference between the range of the continued sums 
of 2 m N tosses of a single coin and N tosses of 2 m coins; 
Xd. = continued sums, or cumulative totals; 
base of natural logarithms; 
number of gains; 
total number of instances in which the losses exceed the gains; 
number of heads r in N throws;h = rN; 
R/e.. 
log N/2’ 
number of losses; 
mean annual draft of water from storage; 
number of coins in any toss NV; »C; (see C); 
number of annual discharges recorded; period of observation; also, 
number of tosses, of a set of m coins; number of departures from 
a mean; 
Nm; 
discharge, annual; 
range, from maximum to minimum, of the curve obtained by plotting 
the cumulative totals of the departures from the mean 2d,, of 
the annual discharges of a stream, taken in order; 
number of heads in m tosses of a coin; also, number of heads; mC; 
(see C); r1,2 = coefficient of correlation between two variables; 
storage required to maintain a draft B less than the mean M; maxi- 
mum storage deficit; 
= number of tails, m — r in m throws or tosses; 
= abscissas; 
= ordinates; 
= “sum of”; and 
= standard deviation: 
a, = o of the number of heads or tails; 


o4 = o of the number of heads minus the number of tails—that 
is, 2r —2m. 
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